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We utilize lattice simulations of the dimensionally reduced effective field theory (EQCD) to de- 
termine the quark number susceptibility of QCD at high temperature {T > 2Tc). We also use 
analytic continuation to obtain results at finite density. The results extrapolate well from known 
perturbative expansion (accurate in extremely high temperatures) to 4d lower temperature lattice 
data. 
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1. Introduction 



The quark number (baryon number) susceptibility is an observable which is one of the sig- 
nature of the quark-gluon plasma in heavy ion collision experiments Thus, it is of interest 
to accurately calculate the quark number susceptibility theoretically. In weak-coupling perturba- 
tion theory, the susceptibility of the quark-gluon plasma has been calculated up to order g^\n\/g 
[0]. Because of the asymptotic freedom, at high enough temperatures the perturbation theory is a 
valid approach. However, the convergence of the perturbative series is bad at physically accessible 
temperatures, and the applicability of the results is not obvious, see Fig. [l[ The order-by-order 
behaviour of the susceptibility is not systematic, and the low-temperature behaviour changes qual- 
itatively at orders and g^ . Further, if we allow variations in the unknown C?(g^) -coefficient in 
the expansion (Fig. [I] right), we obtain quite large variation in the final result at temperatures under 

lor,. 

The quark number susceptibility has also been studied with lattice simulations [^, While 
the standard lattice methods (with dynamical fermions) are the best method to study QCD in the 
immediate vicinity of the phase transition, at higher temperatures it is more economical to simulate 
dimensionally reduced effective theory, electrostatic QCD (EQCD) |j5|]. This method has already 
been used to calculate the pressure in quark-gluon plasma [|6|]. Here we present the updated results 
of simulations of quark number susceptibility for A/j = 2 for zero chemical potential. We also 
present preliminary results from simulations extended to finite chemical potential. These results are 
obtained by doing simulations with imaginary values of the chemical potential and then analytically 
continuing to the real values. This is achieved by fitting a polynomial of /^^ to the data. 
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Figure 1: Left: The perturbative expansion of quark number susceptibility order by order for Nf — 2. The 
coefficient at order has been fixed here to match the lattice measurements. Right: The effect of changing 

the value of the unknown 0{g^) coefficient, parametrised as C{Nf) ■ The perturbative results are 

from [||] and the lattice results from [||, Q]. 
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2. Susceptibility in electrostatic QCD 

EQCD is defined by tlie action 

5e = I d3^|^Tr[i^.2] +Tr[D,-,Ao]2 + m2Tr[A2] + /73Tr[Af,] +A3(Tr[A2])2| , (2.1) 

wliereF,y = diAj — djAi + igj,[Ai,Aj\ andD,- = di + igj,Ai. Fij,Ai and Aq are traceless 3x3 Hermitean 
matrices (Aq = AqT^, etc). The coupling and the mass parameters g3,m^, 73 and A3 are determined 
by the physical 4d temperature, renormalization scale A^, chemical potential jJ. and the number 
of massless fermions. It is convenient to use the dimensionless ratios 

83 83 8l 

which determine the physical properties of EQCD. The -dependence of the parameters is, at 1- 
loop level, 

y = y,-o(^i+y^j^^^, z = lg, x = .,.o, (2.3) 

where p. = jj./ {nT) and the /i = expressions can be found in ref. [|5|]. The two loop corrections 
have been calculated in ref. [^], but the effects remain in practice negligible. 

The quantity we are interested in is the quark number susceptibility, which we define as a 
derivative over one flavor u only: 

Substituting 5e from ( ^ ) we arrive at the result 

Z3 = -2^^.^o(TrAg) 



1 

36 



dh,dh2 ((TrA^(ri)TrA^(r2)) - (TrA^)^) 

2M«3'^=o / d^ndh2 ((TrA2(n)TrA2(r2)) - (TrA^)^) . (2.5) 



V{lNc+Nf) 

Thus, the quark number susceptibility is obtained by measuring the condensates {TrA^), ((TrAQ)^) 
and ((TrAQ)2) on the lattice. Due to the superrenormalizable nature of the theory, measurements 
can be rigorously converted to MS scheme in the lattice continuum limit, and because MS was used 
in the perturbative matching to 4d QCD, this also allows us to compare to 4d results. 

The lattice counterterms needed for the continuum limit of {TrA^) are given in if], and of 
V(TrAQ(ri)TrAQ(r2)) in [^. The contribution including (TrAQ(ri)TrAQ(r2)) is not UV divergent 
and thus does not require counterterms. 

Finally, the relation between and the true 4d susceptibility is 

8% 

X = Yi^3 + J^^P^ (2.6) 

where Ap = /jqcd — P3d is the perturbative 3d^4d matching coefficient for pressure, and can be 
found in ||9|l. 
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Figure 2: Continuum extrapolations of y(((Ti-Ag)2) - (TrAg)^) and V (((TrA^)^) - (TrA^)2) 



3. Lattice measurements 



Lattice simulations are carried out for A'f = 2. We use 6 different values of chemical potential 
/I, and for each value of pL we use eight different values of temperature T . For each of these {pL,T)- 
pairs we use five different values of the lattice spacing a, in order to obtain a reliable continuum 
limit. To check the finite volume effects we did simulations with different volumes at the smallest 
lattice spacing; for a detailed analysis in a related theory see 

Precise continuum limits are necessary for accurate determination of For the condensate 
(TrAg) we use a fit ansatz of form 



+ + 



log(iS) + - 



(3.1) 



j8 j8 j82' 

where j3 = 6/(^30'). The existence of the logarithmic term increases the final errors significantly. 
However, the coefficients of the above ansatz are perturbatively calculable, and there is an ongoing 



program to determine the coefficient of the log-term using stochastic perturbation theory [|1 Ip. The 
knowledge of this would reduce the errors by order of magnitude. 

The contributions V{{{TxAlf) - (TrA^)^) and ^(((TrA^)^) - (TrA^)^) are fitted with second 
order polynomial ansatz 

^1 + ^ + ^- (3.2) 
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Figure 3: Simulations done with two different (x,3')-pairs while varying the imaginary chemical potential 
in the range < —iz < 0.15. Within our accuracy varying iz does change the results, making analytic 
continuation straightforward. 



This fits the data well, see Fig. ^ 

We study the chemical potential dependence by using imaginary jx and performing analytical 
continuation. However, this turns out to be rather trivial: for fixed x,y the dependence of the results 
on iz °^ in is very small and not visible within our statistical errors, see Fig. |3[ (However, see the 
note above about the statistical errors in (Aq).) Thus, the /^-dependence of the results is, in practice, 
completely due to the jtt -dependence of y. 

The final continuum extrapolated results agree well with the perturbative susceptibility. It is 
of the form 

;t:pert = aiy^^^ + a2y + a^y^l'^ + a^. (3.3) 

Hence the difference of lattice and perturbation theory should behave as j^^/^, which is the case, 
as can be seen in Fig. |3[ 
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Figure 4: Left: The quark number susceptibility in EQCD, with different values of chemical potential in di- 
mensionless units. The solid lines are the perturbative result. Right: The difference between the perturbative 
and lattice results. (The T-scale on the top of the figures corresponds to z = case.) 
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Figure 5: The susceptibiUty in 4d units at /i = 0. The results agree with 4d lattice simulation results. 
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After the matching to 4d QCD we obtain the susceptibility in physical units. Our results 
significantly deviate downwards from the perturbation theory, as can be seen from Fig. bringing 
the results closer to the recent simulations by Karsch et al. [Q]. However, one should bear in 
mind that our results suffer from a matching ambiguity related to the unknown 0{g^) coefficient 
in perturbation theory, see Fig. [l], and the results need to be matched to a known point (4d lattice 
simulation) at some low temperature. Nevertheless, we can say that the deviation from perturbative 
result is still rather large at \OTc. See Fig.|[ 
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